Two-photon width and gluonic component of cr//o(600) 
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We analyse data on tt-k and 77 scattering below 700 MeV within an improved analytic K-matrix model. This 
model is based on an effective theory with couplings between resonances, hadrons and photons. The two-photon 
decay of a resonance can proceed through intermediate transition into charged hadrons (here: 7v + tv~) and their 
subsequent annihilation or through a "direct" transition into photons. Our analysis confirms the rather large 
total radiative width of the a resonance which we find as (3.9 ± 0.6) keV but suggests its dominance by the tttv 
rescattering process. This process is not sensitive to the internal structure of the resonance contrary to the direct 
component which we find small, (0.13 ± 0.05) keV, and well consistent with the expectations for an unmixed 
glueball according to the QCD sum rule calculations. 



1. INTRODUCTION 

The study of scalar mesons and the interpre- 
tation of the experimental results is, despite a 
long lasting effort, still an active field of research 
with controversies on experimental results and 
the theoretical interpretation. The lightest scalar 
mesons have been interpreted as conventional qq, 
but also as tetra-quark or molecular states. In ad- 
dition, there is the definitive expectation within 
QCD of the existence of scalar glueballs which 
can mix with the nearby states with quark con- 
stituents. Depending on the mass of the glueballs, 
this affects strongly the interpretation of the spec- 
trum and leads to different scenarios. 

The phenomenological analysis attempts to 
group the spectrum of observed states into appro- 
priate qq or 4q multiplets: left over states are pos- 
sible candidates for glueballs. Also one expects 
the production of such states to be enhanced in 
"gluon rich processes" , while they should be sup- 
pressed in 77 reactions and corresponding rules 
hold for decays as well. 

The existence of glueballs has been predicted 
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long ago in the early time of QCD as consequence 
of the self-interaction of gluons and first scenarios 
have been developed already in 1975 pQ. Today 
there is agreement on the existence of such states 
and the lightest state to be a scalar meson. Quan- 
titative results are available from Lattice QCD 
and from QCD sum rules (QSSR). The theoreti- 
cal results, recalled below, suggest a light scalar 
glueball with mass around 1 GeV or below. 

Among the light particles the /o (600) / a meson 
could be such a gluonic resonance. Recent analy- 
ses of the 77 — > 7T7T processes have extracted the 
width of / (600)/cr -> 77 as: (4.1 ± 0.3) keV g], 
(3.5 ± 0.7) or (2.4 ± 0.5) keV [3] to (1.8 ± 0.4) 
keV [I], while the one from nucleon electromag- 
netic polarizabilities has given (1.2 ±0.4) keV [5]. 
Such results have been interpreted in [5] as dis- 
favouring a gluonic nature which is expected to 
have a small coupling to 77 |6l7l8l9ll0lllj . 

In our recent paper [12] a resolution of this con- 
flict has been suggested. From the analysis of the 
77 — > tttt processes in the low energy region below 
700 MeV we conclude that they are dominated 
by the coupling of the photons to charged pions 
and their rescattering, which therefore can hide 
any direct coupling of the photons to the scalar 
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mesons. Some first results from this study have 
been presented elsewhere |13j . 

2. A LIGHT 0++ GLUEBALL 

There are several results suggesting a light 
scalar glueball: 

• Lattice QCD. Calculations in the simplified 
world without quark pair creation (quenched ap- 
proximation) find the lightest state at a mass 
around 1600 MeV (recent review in Ref. [14]). 
These findings lead to the construction of mod- 
els where the lightest glueball/gluonium mixes 
with other mesons in the nearby mass range of 
around 1300-1800 MeV. However, recent results 
beyond this quenched approximation [15] indicate 
the lightest state with a large gluonic compo- 
nent to fall into the region around 1 GeV, and 
therefore, a scheme based on the mixing of me- 
son states with all masses higher than 1300 MeV 
could be insufficient to represent the gluonic de- 
grees of freedom in the meson spectrum. Fur- 
ther studies concerning the dependence on lattice 
spacings and the quark mass appear important. 

• QCD spectral sum rules (QSSR) and Low- 
energy theorems (LET). These approaches have 
given quantitative estimates for the mass and de- 
cay properties of the gluonic bound states. In 
particular, in a combined analysis of subtracted 
and unsubtracted sum rules a low mass for the 
bare unmixed gluonium state is obtained |6l8j : 

M aB ~ (0.95 - 1.10) GeV , (1) 

besides the heavier state around Mq ~ (1.5 ~ 
1.6) GeV. These masses are similar to the un- 
quenched and quenched lattice results respec- 
tively, quoted above. The hadronic width in this 
approach is obtained rather large: 

T„ B ^ +7! - ~ 0.7 GeV , (2) 

whereas the radiative decay width is found small: 

r (TB ^ 77 ~ (0.2 - 0.6) keV . (3) 

Indeed, this width is considerably smaller than 
the value of about 5 keV, which is expected for a 
state 5*2 = (m + dd) with mass of 1 GeV [TT] . On 
the other hand, QSSR predicts, for a four-quark 



state having the same mass of 1 GeV, a 77 width 
of about 0.4 eV [TT]. 

• Phenomenological studies. A full under- 
standing of the scalar meson spectrum is re- 
quired if the glueball is to be found as a left 
over state after the identification of flavour mul- 
tiplets. Several schemes exist, motivated by the 
quenched lattice result, where the extra gluonic 
state is assumed to mix into the three isoscalars 
/o(1370), / (1500) and / (1710) [IB]. In an al- 
ternative approach |17j . the lightest qq multi- 
ple is formed from / (980), / (1500), 1^(1430) 
and ao(1450), and the glueball is identified as 
the broad object at smaller mass represented by 
/o(600); it dominates irir scattering near 1 GeV 
but extends from threshold up to ~ 1800 MeV 
(see also ref. [18 ). The appearence of this broad 
object in most gluon rich processes was consid- 
ered in support of this hypothesis. In this analysis 
of the spectrum, results from elastic and inelas- 
tic 7T7T scattering as well as from D, B and J/tp 
decays have been considered [17I19) . 

3. THE DIRECT CONTRIBUTION 

The amplitudes for the reactions 77 — > nir 
are largely determined by elastic nn scattering 
and likewise in the case of multi-channel pro- 
cesses. First, there are the extended unitarity 
relations corresponding to the Watson theorem 
for the single channel; secondly, there are dis- 
persion relations where the necessary subtrac- 
tion terms generate polynomial ambiguities 20J. 
This general formalism has been applied by Men- 
nessier [21] . to obtain the electromagnetic pro- 
cesses 77 — > 7T7T, KK, given the strong processes 
7T7T — * 7T7T, KK. The latter ones are represented 
by a K matrix model which represents the ampli- 
tudes by a set of resonance poles. In that case, 
the dispersion relations in the multi-channel case 
can be solved explicitly, which is not possible oth- 
erwise. This model can be reproduced by a set 
of Feynman diagrams, including resonance (bare) 
couplings to 7T7T and KK, and 4-point 7T7t and 
KK interaction vertices. A subclass of bubble 
pion loop diagrams including resonance poles in 
the s-channel are resummed (unitarized Born). 

The radiative width of the resonances cannot 
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be predicted due to the polynomial ambiguity, 
which can be taken into account by introducing 
as free parameters the "direct couplings" of the 
resonances to 77 in an effective interaction ver- 
tex. In the present analysis we have extended 
the model by the introduction of a shape function 
which takes explicitly into account left-handed 
cut singularities for the strong interaction ampli- 
tude. This allows a more flexible parametrisation 
of the 7T7T data at low energies and improves the 
high energy behaviour. The model shows differ- 
ent characteristics at low and high energies: 

• Low energy limit: rescattering. A striking 
feature of the low energy 77 — > irir scattering 
is the dominance of the charged over the neu- 
tral 7T7T cross section by an order of magnitude 
which can be explained by the contribution of 
the pion-exchange Born term in 77 — > it + it~ . In 
the process 77 — > ir ir , the photons cannot cou- 
ple "directly" to 7r°7r° but through intermediate 
charged pions and subsequent rescattering with 
charge exchange. This feature is realized in the 
analytic model considered here [21] which has it 
in common with the Chiral perturbation theory 

• High energy limit: direct production. On the 
other hand, at high energies the photon can re- 
solve the constituents of the hadrons. An example 
is the production of /2(1270) in 77 — > tttt. The 
analytic model [21] predicts a decreasing cross 
section for this process with pion exchange reach- 
ing < 10% of the observed cross section at the 
peak of the fi . These considerations require the 
need of a "direct coupling" of fi — > 77. Indeed, it 
is well known that the radiative decays of the ten- 
sor mesons fi , , a,i are well described by a model 
with direct coupling to the qq constituents ac- 
cording to the SU(3) structure of the nonet with 
nearly ideal mixing (see e.g. Ref. [23]). We there- 
fore interpret the direct terms in the 77 processes 
as originating from the coupling of the photon to 
the partonic constituents of the resonances (such 
as qq, 4q, gg,. . .). 

4. ANALYTIC K-MATRIX MODEL 

In the application of this model [21] . we first 
obtain a suitable K-matrix parametrisation of 



the 7T7T scattering data, and then determine the 
direct coupling by comparing the model with 
the 77 results. In the present analysis, we re- 
strict ourselves to the low mass region below 700 
MeV where we neglect vector and axial-vector 
exchanges [21], inelastic channels and D- waves. 
Furthermore, we assume a pointlike pion-photon 
coupling, which is expected to be a good approx- 
imation in the region where we are working. 

We apply the analytic model [21J, but we intro- 
duce a shape function fo(s) which multiplies the 
(77771" coupling. The real analytic function fo(s) is 
regular for s > and has a left cut for s < 0. For 
our low energy approach, a convenient approxi- 
mation, which allows for a zero at s = sa and a 
pole at an > simulating the left hand cut, is: 



Ma) = 



s - SAO 



(4) 



s + o- D0 

For simplicity, we don't include the 4-point cou- 
pling term. The unitary 7T7t amplitude for 
the isospin 1 = S-wave is written as: 

Gfo(s) _ Gfo(s) 



SR 



Gf(s) -D(s) 



(5) 



where G = g^ B is the bare coupling squared. 
Unitarity determines the imaginary part, while 
a dispersion relation subtracted at s = deter- 
mines the form of f(s) [12\ . 

Using the S wave amplitude in Eq. ([5]) we de- 
rive the amplitude Ty for the electromagnetic 
process for isospin / = as: 

if +<*#)+«!■ ro 

The first contribution comes from the Born term 
for 77 — * 7r + 7r~, the second one from the 7T7t 
rescattering and the third one represents the poly- 
nomial contribution V = sF 1 \[2 which reflects 
the ambiguity from the dispersion relations and 
represents the direct coupling of the resonance to 
77. The residues at the pole sq of the rescattering 
and direct contributions to in Eq. © deter- 
mine the respective branching ratios. Similarly, 

(2) 

we construct the 1 = 2 5-wave amplitude T 

(2) 

and as well as the cross sections for the 7T7t 
and 77 scattering processes [2"T] . 
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Figure 1. 7r°7r° cross-section [nb] versus \fs [GeV] 
and fit by analytical model [21]: F 7 = (dot- 
dashed); F 7 = -0.09: 1=0 (large dashed), 1=0 
+2 (continuous); F 7 = -0.07: 1=0+2 (small da- 
shed). Data from Crystal Ball [27J (| cos 6»| < 0.8). 
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Figure 2. Same as in Fig. [T]but for ir + ir : F 1 = 
(dot-dashed); F 1 = -0.07 and 1=0+2 (small 
dashed). The continuous line corresponds to the 
non-unitarizcd Born amplitude with F 7 = 0. 
Data are MARKII [2J (| cos6»| < 0.6). 

5. RESULTS OF THE ANALYSIS 

We first determine the parameters of the model 
for 7T7T scattering below 700 MeV from the best 
approximation of our formula to the phase shifts 
<5g°' ) and 5^ obtained in the Roy equation analy- 
sis in [24]. For both channels 1 = and / = 2 we 
fit 4 parameters each. We find the pole mass: 

M a ~ (422 - i 290) MeV (7) 

which is close to the masses in other recent deter- 
minations (441 - i 272) MeV [25] or (489 - i 264) 
MeV [25]. 



Given the tttt amplitude, we can predict the 
cross sections for 77 — > tttt where the only free 
parameter F 1 is related to the strength of the 
direct coupling a — > 77. The fit of the model 
to the data from Crystal Ball (7r°7r°) [27 and 
MARK-II (7r + 7T~) [28] collaborations is shown in 
Figs. [T] and [2] for the cases with direct contribu- 
tion and without ("unitarised Born" : F y = 0). 
We obtained excellent fits to the neutral pion data 
with F 7 = —0.090. A fit to both channels yields 
F 7 = —0.070. Then the fit in the charged chan- 
nel deviates from the data at high-mass; the large 
systematic errors and the absence of data points 
in the region 0.40 to 0.55 GeV do not permit a 
good understanding of this channel. We take as 
a final estimate F 1 ~ -(0.080 + 0.014). From the 
residues at the pole sq and r^—,^ — 2ImM cr , we 
can deduce the "partial" 77 widths at the com- 
plex pole : 

r^ 77 (0.13 + 0.05) keV , 
r^ c 77 ~ (2.7+0.4) keV , (8) 
and the total 77 width (direct + rescattering): 
~ (3.9 ± 0.6) keV . (9) 

Our result for the direct coupling is comparable 
with our previous value [13] obtained without us- 
ing the shape function fo(s), but not as small as 
in a parallel analysis [29] , Our total 77 width is 
compatible with the range of values obtained pre- 
viously [2l3l4l5j . These analyses, however, don't 
separate the direct term. 

The results from QSSR/LET are obtained in 
the physical region and can better be compared 
with experiment by using the corresponding re- 
sults for the "visible meson" on the real axis in- 
stead of the results at the complex pole: this 
can be, either the Breit-Wigner mass and width 
(see 13 ) or the "on shell" mass and width (see 
e.g. [30J), which are obtained for the mass M° s 
where the real part of the propagator vanishes 
ReX»((M° s ) 2 ) = 0. For this definition, we obtain: 

M° s w 0.92 GeV, ~ !-02 GeV, 

m (1.0 + 0.4) keV. (10) 

As the model is extrapolated here towards ener- 
gies beyond its validity (M ~ 1 GeV), we consider 
these results as a crude approximation. 
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6. COMPARISON WITH QSSR/LET 

One can notice that the QSSR and LET pre- 
dictions for the mass and hadronic widths of a 
low mass gluonium ctb in Eqs. ([T]) and are 
in remarkable agreement with the results in Eq. 
(fTOf for an on-shell resonance. The direct 77 cou- 
pling, which can reveal the photon coupling to 
the intrinsic quark (qq, Aq, . . .) or gluon ( gg, . . .) 
constituent structure of the resonance, can be re- 
lated to the QSSR/LET evaluations of its width 
through quark or gluon loops as in case of the 
quark triangle for the pion or the /a(1270). Our 
a — ► 77 width in Eq. (fit))) (also the pole value 
in Eq. ([5])) are consistent with the small value 
predicted in Eq. ([3]). This "overall agreement" 
favours a large gluon component in the c//o(600) 
wave function. It disfavours, in particular, a 
qq interpretation because of the larger radiative 
width to be expected (~ 5 keV). On the other 
hand, a four-quark interpretation would require 
a smaller radiative width and is disfavoured as 
well by our result in Eq. (at the level of 2.5 
a) . 

Improvements of our estimates need more pre- 
cise data below 700 MeV, which could be pro- 
vided by the KLOE-2 experiment in the future 
31j. Furthermore, an extension of the analysis to 
higher energies will be important. 
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